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PROGRAM STRUCTURE 


Swot is constructed and formulated in such a manner that the sequential 
flow of program control necessary for solution using high speed digital com- 
puters can be easily understood. The basic computer language used is Fortran V 
and the model has been successfully applied using the CDC 6600 and 6400, the 
UNIVAC 1108 and 1106 and the RCA Spectra 70/45 computers. In its present form 
HYDTID is essentially machine independent. The computation time and storage 
depend on the size of the system being modeled, the mesh size and time step 
being used, the number of computational or water cells in the grid network, and 
the length of simulation time desired. Recent applications of the model to 
Texas bays have required as much as 150,000 words 0 and from 2 to 6 
minutes of UNIVAC 1108 computer time to simulate one 25 hour tidal cycl@e> For 
the Masonboro Inlet problem, the coarse grid model required about four minutes 
of UNIVAC 1108 time to simulate one 12.5 hour tidal cycle and to generate the 
input flows for the fine grid sub-model. The fine grid sub-model required about 


40 minutes of computer time to simulate one cycle. — 


Basically the model consists of an executive control program and eight 
subroutine packages for performing specific computational tasks. Input data 


to the program are read from cards or magnetic tape, and output is either 


printed, punched on cards, or written on magnetic tape. Detailed descriptions 


of the basic program elements are presented in the following sections. 


Executive Control Program. This is the basic command element of Program 
HYDTID. All program control data, real system description data, and model 


operation data required for computing tidal hydrodynamics are input into the 
program by this routine. Variables are initialized, constants are set, and 

a portion of the basic input data is printed out. At every time step, control 
is transferred from this routine to appropriate subroutines for calculation of 
instantaneous values of tidal amplitudes and flows, as well as net flows, 

net velocities and mean velocities. Input data are processed appropriately 
based on whether the program is being applied to the fine or coarse grid con- 
figuration. This element of HYDTID also transfers control to appropriate 
subroutines for printing tidal amplitudes and discharges for selected cells 

at specified time intervals, for storing instantaneous hydrodynamics for all 
cells on magnetic tape at specified time intervals, for storing tidal amplitudes 
for selected cells for plotting, and for storing final values of the basic 


variables on cards or tape at the end of program execution. 


Subroutine CALTID. This subroutine computes at every time step the 
appropriate values of the four exciting tides that are specified around the 
periphery of the coarse grid model to impart tidal fluctuations in the system. 
Hourly values of tidal amplitudes read into the program from cards are interpo- = 


lated appropriately to ootain tidal values at every time step. 


Subroutine PRINTI. This subroutine prints the remaining basic input 
data not printed in the Executive Control Routine. Appropriate descriptive 
headings and titles are printed with the data so that it is possible to check 
that all prototype conditions are properly accounted for in the model. Print 
out of program control parameters serves as a check to see that the program 
has operated correctly. Proper specification of certain program control 
parameters can cause control to pass over this subroutine or execute only a 


portion of it. 


Subroutine CALCQH. This subroutine is the basic computational element 
of the tidal hydrodynamics model. Control must pass through this subroutine 
during each time step of the computation process. In this subroutine, control 
moves from cell to cell according to a previously defined sequencing scheme so 
that only water cells are considered in the calculations. Based on the type 
of boundary condition required by a particular computational grid cell as 
specified by its identifying "flag" value, control is transferred to appropriate 
statements within this subroutine for calculation of certain constants and 
coefficients. Utilizing these constants and coefficients and known values of 
tidal amplitudes and flows from the previous time step, new values of these 
quantities are then determined explicitly at the end of the routine. One 
execution of this subroutine provides a complete new array of tidal amplitudes 
and flows per foot of width in the two coordinate directions. Water depths 


are computed from the new values of tidal amplitude at the end of this routine. 


Subroutine PRINTO. This subroutine outputs the results of the basic 
model computations at specified time intervals. Specifically, values of tidal 
amplitudes and discharges are printed out for twenty pre-specified grid cells 
located in the grid system. Ordinally, hourly values are printed out, 
however, any desired time interval greater than the computational time step 
can be used. In the event that the complete array of final computed values 
of tidal amplitudes and flows per foot of width for all grid elements are 
required to be saved at the end of model operation, control can be transfer- 
red to this subroutine where these values are either punched on cards or 
written on magnetic tape. These ending values of the basic variables are 
used as initial values for subsequent computer runs. Also at the end of model 
operation, hourly velocities for the above 20 cells are punched for the x- and 


y-directions. 


Subroutine NETVQD. This is an optional subroutine that calculates net 
velocities, net flows, or average depths which occur during a tidal cycle for 
all computational grid elements. Net velocities and flows and average depths 
are printed for all grid cells and can also be punched on cards or stored on 


tape by this routine. 


Subroutine STRVEL. This is also an optional subroutine, and it stores 
on magnetic tape the instantaneous hydrodynamics for all cells at specified 
time intervals. This routine is particularly useful for storing hydrodynamics 
at times of ebb, flood, and slack tides so that flow conditions throughout the 
system can be analyzed at these times. Another use of this routine has been 
to store instantaneous hydrodynamics for all cells at short time intervals 


and to use these data in a marker particle program to trace paths traversed 


by water particles during a tidal period. 


Subroutine PLOTHS. This subroutine is an optional plot package which 
when executed, results in verification type plots of tidal amplitude for 
specified grid cells. Both the tidal amplitude computed by the hydrodynamic 
model and the tidal amplitude measured in the prototype are plotted for 
specific locations. Appropriate descriptive labels and titles are also 
included on the plots. As many as twenty different locations (grid cells) 
can be chosen for tidal amplitude plots. By obtaining these plots for several 
locations throughout the system being simulated, an idea of the accuracy of 


the model can be obtained. 


Subroutine RITAP. This subroutine is used only in the operation of 
the coarse grid model and compiles arrays of selected flows and tides from 


the coarse grid model results. These selected flows and tides are then 


interpolated temporally and distributed spatially to obtain the boundary input 
flows for the fine grid sub-model. The fine grid flows are then stored on 
magnetic tape for use in operating the fine grid sub-model. For the Masonboro 
Inlet problem, there are two versions of Subroutine RITAP which must be inter- 
changed in the HYDTID program deck depending on the particular inlet geometry 
being simulated. Version I applies to the pre-project condition without the 
jetty, and Version II is used for the post-project condition with the jetty 

in place. Two different versions are needed because the configuration of 
boundary flow cells used to excite the fine grid sub-model are different for 
the two inlet conditions, and therefore different coarse grid flows must be 


interpolated and distributed accordingly. 


DEFINITION OF PROGRAM VARIABLES 


Program HYDTIO 


ANGCOR 
CB 
CELSID 
CODE 


CON] 
CON2 
CT 


DT 
DTODS 
DTO2DS 
OTS 


- Angle between north and x-axis measured clockwise from north. 


- Submerged or overtopping barrier coefficient. 


Literal description equal to SIDE or TOP. 


- Literal, (CARD, TAPE, NONE, BOTH) which designates mode 
of I/0. 


- Base value counter for tidal curve interpolation. 

- Base value counter for tidal curve interpolation. 

- Tidal discharge coefficient. 

- Total water depth in a given cell. 

- Genera] purpose input variable used for temporary data storage. 
- Cell side dimension. 

- Computational time step. 

- DT/DS. 

DT/(2-DS). 


- DT/2. 

- General purpose input variable used for temporary data storage. 
- Rate of evaporation. 

- Literal which denotes end of input data file. 

- Literal which denotes end of input title file. 

- Manning's "n" bottom roughness coefficient. 


- Function of Manning's "n" for computations in the x-direction. 


RY 

G 

GC 
GCDT04 
GDTODS 
GTIDE 


Gl 
G41 


HN 
HPLT 
HPRT 
HPRTA 
HSHIFT 
I 

IBAR 
IBASIC 


ICLL 


IDCARD 
IOTIDE 


Function of Manning's "n" for computations in the y-direction. 


Acceleration of gravity. 

Internal computation constant. 
Internal computation constant. 
Internal computation constant. 


Current value of exciting tide for one of four input tidal 
conditions. 


Exciting tidal elevation temporarily stored for printed output. 
Exciting tidal elevation temporarily stored for printed output. 
Exciting tidal elevation temporarily stored for printed output. 
Exciting tidal elevation temporarily stored for printed output. 
Current tidal elevation in a given cell. 

Prototype tidal elevation used in verification plots. 

Newly computed tidal elevation in a given cell. 

Storage variable equal to tidal elevation to be plotted. 


Print out variable for tidal elevation. 


Print out variable for tidal elevation. 

Elevation difference between MSL and datum of input data. 
Standard grid column indicator. 

Grid column indicator for submerged barrier cells. 


Internal variable which indicates number of program options 
desired. 


Grid column indicator for water cells where computations are 
required. 


Variable which indicates mode of basic cell data input. 


Identification number which assigns exciting tide to appropriate 
cells. 


IDUM General purpose input variable used for temporary data storage. 
IFLAG Computational cell flag number which denotes type of calculations 

to be performed. | 
IFLOW Grid column indicator for external inflow cells. 
IHKP Grid column indicator for cells in coarse grid where ending 

H-values are to be punched for input into fine grid sub-model. | 
IK Internal counter. 
ILB Internal counter. 
ILF Internal counter. 
IMAX Total number of columns in grid. 
IMXJMX Total number of cells in grid. 
INETFL Variable which specifies net flow option. 
INEW Internal variable used to facilitate I/0. 
IODISP Variable which specifies mode of dispersion coefficient output. 
TONFLO Variable which specifies mode of net flow output. 
IONVEL Variable which specifies mode of net velocity output. 
IP Grid column indicator for cells where tides and flows are to 

be periodically printed. 
IPDATA Variable which denotes extent of input data print out. 
IQHIN Variable which specifies mode of initial hydrodynamics input. : 
ISAVQH Variable which specifies model of final hydrodynamics output. 
ITIDE Grid column indicator for tidal excitation cells. | 
IVLTAP Variable which specifies mode of instantaneous hydrodynamics 

output. 
J Standard grid row indicator. 
JBAR Grid row indicator for submerged barrier cells. 
JCLL Grid row indicator for water cells where computations are required. 
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JFLAG 


JFLOW 
JHKP 


KD 
KEPSAV 


KG 
KINDAT 
KINIQH 
KK 
KQ 


KODISP 
KONETF 
KONETV 


KOTVEL 


Two digit cell flag which specifies the particular finite 
difference formulation of the convective acceleration cross- 
product term that is to be used for a given cell. 

Grid row indicator for external inflow cells. 


Grid row indicator for cells in coarse grid where ending 
H-values are to be punched for input into fine grid sub-model. 


Internal counter. 
Internal counter. 


Internal counter. 


Total number of rows in grid. 


Grid row indicator for cells where tides and flows are to 
be periodically printed. 


Grid row indicator for tidal excitation cells. 


Internal counter. 


Temporary counter for submerged barriers. 


Temporary counter for external inflows. 


Temporary variable used to indicate storage of hydrodynamics at 


end of one tidal 


Internal counter. 


Tape unit number 
Tape unit number 
Internal counter 


Internal counter 
selected cells. 


Tape unit number 
Tape unit number 


Tape unit number 


Tape unit number 


cycle. 


for reading basic cell data. 
for reading initial hydrodynamics. 
used in data input. 


used in printing basic hydrodynamics for 


used for storing dispersion coefficients. 
used for storing net flows. 
used for storing net velocities. 


used for storing instantaneous hydrodynamics. 


KOUNT 


KOUTDA 
KPRINT 


KQCTP 


KQFTP 


KRSOFN 
KT 


LINMAX 


M 

MA 

N 
NFLOW 
NN 
NPLOT 
NPRPLT 


NREEF 
NTIDE 
OMEGA 


PI 


Internal counter used to designate specific water cells. 
Tape unit number for storing final hydrodynamics. 


Variable which controls punching of hydrodynamics at end of 
one tidal cycle. 


Tape unit number used for storing selected flows from coarse 
grid model for input to fine grid sub-model. 


Tape unit number used for storing external inflows for exciting 
fine grid sub-model. 


Variable which indicates type of model operation to be performed. 
Internal counter. 


Variable which indicates number of sets of hydrodynamic output 
to be punched per page. 


Internal counter for plotting tides. 

Internal counter. 

Internal counter. 

Total number of external inflows. 

Internal counter. 

Total number of cells where tidal plots are to be made. 


Variable which designates the order of 20 specified cells where 
basic hydrodynamics are to be periodically printed. 


Total number of submerged and overtopping barriers. 

Total number of external tidal excitation cells. 

Coriolis parameter. 

Constant equal to m (3.1416). 

Time interval for printing basic hydrodynamics at selected cells. 
External inflow for a given cell. 


Current value of flow per foot of width in x-direction for a 
given cell. 


QXN 


QYN 


REMARK 
SIDE 
SQTG 
STATON 


TCOUNT 
THETA 
THETA] 
TIDE) 
TIDE2 
TIDE3 
TIDE4 
TID1 
TID2 
TID3 
TID4 
TIM 
TIME 
TIMEIN 
TIMTOT 


Newly computed value of flow per foot of width in x-direction 
for a given cell. 


Current value of flow per foot of width in y-direction.for a 
given celi. 


Newly computed value of flow per foot of width in y-direction 
for a given cell. 


Rainfall] rate. 

Variable used for storing title inputs. 

Literal used to designate right side of computational cell. 
Square root of G. 


Literal used to identify specific cells where basic hydro- 
dynamics are periodically printed. 


Time counter used for printing basic hydrodynamics. 

Wind angle. 

Wind angle temporarily stored for print out. 

Tidal elevation read into program for Exciting Tide No. 1. 
Tidal elevation read jnto program for Exciting Tide No. 2. 
Tidal elevation read into program for Exciting Tide No. 3. 
Tidal elevation read into program for Exciting Tide No. 4. 
Current interpolated tidal elevation from Exciting Tide No. 1. 
Current interpolated tidal elevation from Exciting Tide No. 2. 
Current interpolated tidal elevation from Exciting Tide No. 3. 
Current interpolated tidal elevation from Exciting Tide No. 4. 
Time stored for plotting. 

Current value of simulated time during model operation. 
Beginning time of model operation. 


Total time to be simulatec. 


TIMVEL 


TMARK 
TMAX 
TNET 


TOP 
TPER 
TPLOT 


UAPRT 


UAPRTA 


UPLT 


VAPRT 


VAPRTA 


VPLT 


ZB 


Specified time interval for storage of instantaneous hydro- 
dynamics. 


Internal time counter for storage of instantaneous hydrodynamics. 
Final value of time at end of model operation. 


Value of time at wnich computations for net flows and velocities 
begin. 


Literal used to designate top side of computational cell. 
Period of tidal cycle. 


Value of time at which storage of hourly H-values begins for 
tidal plots. 


Variable used for printing flows in x-direction at specified 
cells. 


Variable used for printing flows in y-direction at specified 
cells. 


Velocity in x-direction punched at PTIME intervals for selected 
cells where velocity comparisons are desired. 


Variable used for printing flows in y-direction at specified 
cells. 


Variable used for printing flows in y-direction at specified 
cells. 


Velocity in y-direction punched at PTIME intervals for selected 
cells where velocity comparisons are desired. 


Wind velocity. 

Temporary storage variable generally set equal to W. 
Effective wind stress term for x-direction. 
Effective wind stress term for y-direction. 

MSL elevation of bottom of cell. 


MSL elevation of crest of submerged or overtopping barrier. 


Subroutine 


DELT1 


DELT2 


DELT3 


DELT4 


NTID 
NTIOP1 


CALTID 
- Current incremental change during one time step for Exciting 
Tide No. 1. 


- Current incremental change during one time step for Exciting 
Tide No. 2. 


- Current incremental change during one time step for Exciting 
Tide No. 3. 


- Current incremental change during one time step for Exciting 
Tide No. 4. 


- Counter used in interpolation. 


- Counter used in interpolation. 


Subroutine CALCQH 


COEFX 


COEFY 


DBARX 


DBARY 


DBX 
DBY 
DCON 
HMAX 
IFL 
| IFLG 
JAFL 


- Computed coefficient used in basic hydrodynamic computations 
of flow in x-direction. 


- Computed coefficient used in basic hydrodynamic computations 


of flow in y-direction. 


- Internally computed variable involving depths in adjacent 
cells in x-direction. 


- Internally computed variable invoiving depths in adjacent 
cells in y-direction. 


- Average water depth over submerged barrier on side of cell. 
- Average water depth over submerged barrier on top of cell. 
- Reciprocal of average of depths in adjacent cells. 

- The greater of two adjacent MSL water surface elevations. 

- Temporary variable equal to IFLAG. 

- Temporary variable equal to IFLAG. 


- First digit of JFLAG value. 
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JBFL 
JFL 


KBT 
KTT 
QBARX 


QBARY 


QDIFXS 
QDIFYS 
QXBAR 
QYBAR 
SIGN 


ZMAX 


Second digit of JFLAG value. 


Temporary variable used to indicate type of flow calculations 
required at a particular cell. 


Temporary counter for submerged and overtopping barriers. 
Temporary counter for tidal excitation cells. 


Magnitude of actual velocity vector used in calculation of flows 
in x-direction. 


Magnitude of actual velocity vector used in calculation of flows 
in y-direction. 


Flow gradient in y-direction used to approximate 9q,/ ay. 
Flow gradient in x-direction used to approximate 9q,/2x. 
Average flow in x-direction defined at same location as dy: 
Average flow in y-direction defined at same location as qy° 


Temporary algebraic sign variable which indicates flow direction 
across overtopping barriers. 


The greater of two adjacent cell MSL bottom elevations. 


Subroutine NETVQD 


DEPTH 
DXA 
DYA 
QNETX 
QNETY 
UAVE 


VAVE 


Average water depth in a given cell over a tidal cycle. 

Average of water depths in two adjacent cells in x-direction. 
Average of water depths in two adjacent cells in y-direction. 
Net flow for a given cell in the x-direction over a tidal cycle. 
Net flow for a given cell in the y-direction over a tidal cycle. 


Mean tidal velocity in x-direction for a given cell during one 
tidal cycle. 


Mean tidal velocity in y-direction for a given cell during one 
tidal cycle. 
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DX - Dispersion coefficient in x-direction computed using Random 
Process Analogy. 


DY - Dispersion coefficient in y-direction computed using Random 
Process Analogy. 


VNETX - Net velocity in the x-direction for a given cell during one 
tidal cycle. 
VNETY - Net velocity in the y-direction for a given cell during one 


tidal cycle. 


Subroutine STRVEL 


TAPTIM - Current value of time written on tape for checking purposes. 


Subroutine PLOTHS 


A - A processor variable for plotting. 
ACOLMN - Storage vehicle for a plot character. 
ADOT - Print character DO 

AEQUAL - Print character "=". 

Al - Print character "I". 


AMINUS - Print character "-". 


APLUS - Print character "+". 

ASTRSK - Print character "*”. 

BLANK - Print character "*". 

co - Print character "0". 

DIFHF - Internal processor variable. 
DIFHP - Internal processor variable. 
HF - Prototype tidal elevation. 
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Icc - Internal processor variable. 

IHF - Temporary integer storage for prototype tidal elevation. 

IHPLT - Temporary integer storage for model tidal elevation. 

ITCONT - Internal processor variable. 

ITID - Internal processor variable. 

ITIDM] - Internal processor variable. 

ITIDPR - Internal processor variable. 

MM1 - Internal processor variable. 

TIDPRT - Internal processor variable. 

TITEL - Specified literal title of plot. 

TITELY - Literal ordinate label. 

Subroutine RITAP 

DTOT - Sum of k water depths-used as a proportioning base to distribute 
one coarse grid flow to k fine grid cells. 

HOLD - Value of water elevation for beginning of coarse grid time step. 

HTP - Value of water surface elevation at end of coarse grid time step. 

HT PU - Water surface elevation at intermediate time level used to deter- 
mine fine grid input flows. 

KCT - Internal interpolating counter. 

KCTM - Interpolation factor equal to (coarse grid time step/fine grid 
time step). 

Q - Interpolated and distributed value of external inflow for fine 
grid sub-model. 

QOLD - Coarse grid flow at beginning of coarse grid time step. 

Qs - Computed fine grid external inflow in x-direction. 

QT - Computed fine grid external inflow in y-direction. 
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QTPU 
TIME 
TMAX 
ZT 


Coarse grid flow at end of coarse grid time step. 
- Coarse grid flow at intermediate time level. 

- Internal time counter (seconds). 

- Total real time of model operation (seconds). 


- Input variable of cell bottom elevations. 


DATA INPUT 


The data input structure for HYDTID is dependent on the mode of model 
operation and the various program options the user wishes to employ. For 
purposes of this study, three different types of operation modes are defined 
as follows: (1) Coarse Grid Production Run meaning operation of the coarse 
grid model for the purpose of generating the input flows to the fine grid 
sub-model; (2) Fine Grid Production Run meaning any operation of the fine 
grid sub-model; and (3) Coarse Grid Non-Production Run meaning operation 
of the coarse grid model for purposes other than to generate fine grid sub- 
model inputs. Input data is read from both cards and magnetic tape, with 
some data specified in the program itself. In all there are eleven different 
card data files which can be read, however only six of these are necessary 
for coarse grid model operation, and five are required by the fine grid sub- 
model. The contents of the eleven files are described in the subsequent para- 


graphs followed by their appropriate format structures. 


Title File - Titles for First Page of Output 


Four separate 68 character titles can be specified using this file. 
They appear on the first page of the edited output and can be used to describe 
the various conditions under which the model is being operated. The entire 


Title File is also echo printed at the beginning of each run. 


a ee eo ee so cnn nc mimendaiiessaeasmanameeeenmmmeammmmeatl 


Data File A - Program Control Parameters 


The I/0 mode for various types of data in the model are specified in this 
file by assigning the appropriate literal, CARD, TAPE, BOTH, or NONE, in the 
proper space on the File A cards. If TAPE or BOTH are assigned, the tape unit 
number must also be specified. Based on this information, HYDTID performs the 


necessary I/0 operations. 


Data File B - Basic Model Operation Parameters 


Included in the file are the basic parameters which are used in the model. 
All of the parameters are read as floating point variables from columns 74 
through 80 and then assigned to appropriate variable names in the program. The 
various parameters required are described on the format forms which follow and 


the required units are also spec*fied. 


Data File C - Basic Cell Data 


One data card for grid cell included in a model is read by the program 
from this file. Each card is identified with I and J coordinates and includes 
all of the descriptive data necessary for hydrodynamics to be determined for 
every cel] in the computational grid. These data include the following: 


1.  IFLAG - Computational cell flag determined from the individual 
boundary conditions at the cell . 


Be) - Average bottom or ground elevation (feet) referred to 
same datum specified in Data File B, Card 20. 
oe «OF - Manning's "n* value. 
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4. I[DTIDE 
5.  QINFLO 
6. CBX 
7.  ZBX 
&. BY 
9 ZBY 
10.  NPRPLT 
11. STATON 


Tidal identification number (1 < IDTIDE < 4) which 
assigns appropriate exciting tide to the cell if it 
is flagged accorcingly. Otherwise IDTIDE is zero. 


External inflow magnitude (cfs) if cell is flagged 
accordingly. Otherwise QI\FLO set equal to zero. 

jan must be specified to be consistent with coordi- 
nate axes. For fine grid sub-model, this quantity 
does not need to be specified since external flows are 
read from tape for each of the exciting flow cells. 


Discharge coefficient assigned to barriers parallel 
to x-axis when cell flagged accordingly. Otherwise 
CBX set equal to zero. 


Crest elevation assigned to barriers parallel to 

x-axis when cell flagged accordingly. Otherwise 

ZBX set equal to zero. Referred to same datum specified 
in Data File B, Card 20. 


Discharge coefficient assigned to barriers parallel to 
y-axis when cell flagged accordingly. Otherwise CBY 
set equal to zero. 


Crest elevation assigned to barriers parallel to y- 
axis when cell flagged accordingly. Otherwise ZBY set 
equal to zero. Referred to same datum specified in 
Data File B, Card 20. 


Print/Plot order number assigned to 20 selected cells 
for periodic output of basic hydrodynamics and for 
plotting tidal elevations. Otherwise NPRPLT set equal 
to zero. 


Literal station name used as heading when printing 
basic hydrodynamics for 20 selected cells. Otherwise 
STATON left blank. 


Data File D - Exciting Tides 


As the coarse and fine grid models are currently structured, this data 


file is only required by the coarse grid model. For each of the four exciting 


tides used in the model, 26 hourly values of tides (2 tidal cycles) are read 


preceded by an appropriate title card. 


in Data File B, 


SO LOL LLL LE LLL LO TT 


Card 20. 
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Datum for the tides is the same specified 


Data File E - Cell Identification for Storing H-Values 


This is an optional data file used only for production runs with the 
coarse grid model. I and J coordinates are read for those cells in the coarse 
grid model where tidal elevations at the end of one tidal cycle are required 
to establish the initial water levels in the fine grid sub-model. As currently 


structured the program reads 32 sets of coordinates. 


Data File F - Two-Digit Convective Acceleration Cell Flags 


This data file is required for all modes of model operation. For those 
cells that require finite difference approximations of the flow gradients in 
the cross-product terms that are different from the normal centered difference 
formulations, I and J coordinates and an appropriate two-digit flag are read 
from a single card. When a blank card is encountered, the program assumes 


that the end of this file has been reached and control is transferred ac- 


cordingly. 
Data File G - Complete Array of Initial Hydrodynamics ‘% 


Once the models have been operated for a complete tidal cycle, the 
ending values of the hydrodynamics for every cell can be saved and used as 
the initial conditions at which to begin subsequent simulations with the 
models. In this manner, simulations can be made with several short one 
tidal cycle runs rather than a single long run which might take several 
hours of computer time. These ending hydrodynamics can be obtained on cards 


or magnetic tape by specifying the appropriate option in Data File A for the 
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initial run, and can then be read into the models for subsequent runs as 
Data File G again using the option in Data File A. In the event no initial 


conditions are available, initial hydrodynamics are set equal to zero. 


Data File H - Initial H-Values 


This data file is used only with the fine grid sub-model when initial 
hydrodynamics for all cells are not available from previous runs. For this 
situation, all flows are set equal to zero, but initial water levels are 
established from those computed using the coarse grid model. Data File H 
is punched when the coarse grid model is operated under Mode 1, Production 


Run. 


Data File I - Selected Bottom Elevations for Fine Grid Sub-Model Inflow Cells 


This data file is required only for the operation of the coarse grid 
model under Mode 1, Production Run. Data File I consists of the bottom 
elevations at those fine grid cells where exciting inflows are specified. 
These data are used in spatially distributing the computed flows from the 
coarse grid model to obtain the fine grid inputs. The datum for these ele- 


vations is the same as that specified in Data File B, Card 20. 


Data File J - Prototype HF-Values for Tidal Plots 


When tidal verification plots are desired using either model, the 
measured prototype tidal elevations must be read into the program as Data 
File J. The datum for these elevations is the same as that specified in 


Data File B, Card 20. 
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FLOW CHART FOR MAIN HYDTID PROGRAM 


i 
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START 


} 
SET CONSTANTS 


V) 


READ IN 
PROGRAM 
CONTROL 
DATA 
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Legend for Model Type 
CNP - Coarse Grid Model, 
Non-Production Run 


CP - Coarse Grid Model, 
Production Run 


FP - Fine Grid Model 
Production Run 


READ IN BASIC 
MODEL OPERATION 
PARAMETERS 


READ IN BASIC 
CELL INPUT 
DATA 


READ TIDAL 
EXCITATION DATA 


READ CELLS WHERE 
ENDING H VALUES 
ARE TO BE STORED 
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READ TWO-DIGIT FLAG 
ASSIGNMENTS FOR 
CONVECTIVE ACCEL. 


PRINT PROGRAM OPERATION 
INFORMATION AND MODEL 
CHARACTERISTICS 


PRINT ALL SUBROUTINE 


INPUT DATA PRINTI 
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ESTABLISH 
INITIAL HYDRODYNAMIC 
H, QX, QY 


SET EQUAL TO ZERO 


SOURCE OF 
INITIAL HYDRO- 
DYNAMICS 


Card 


READ INITIAL H VALUES 
DERIVED FROM COARSE 
GRID MODEL 
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INITIALIZE 
VARIABLES 


STORE INITIAL 
HPLT VALUES 
ENTRY HPLOT 
SUBROUTINE PLOTHS 


TIDAL PLOTS 
DESIRED 


STORE INITIAL 
HYDRODYNAMICS 
SUBROUTINE STRVEL 


STORE 
INSTANTANEOUS 
HYDRODYN. 


INCREMENT TIME 
c= tT ae 
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STORE SELECTED FLOWS 
& TIDES FOR DETER- 
MINING FINE GRID INPUTS 


ah 


SUBROUTINE RITAP 
ENTRY RITCTP 
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PRINT HYDRODYNAMICS 
FOR 10 SELECTED CELLS 
if TIME = PRINT TIME 


SUBROUTINE 
PRINTO 


STORE HPLT VALUES 
SUBROUTINE PLOTHS 
ENTRY HPLOT 


OPTIONS 
REQUIRED OTHER 
THAN BASIC 
HYDRODYNAMICS 
COMPUTATION 
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COMPUTE 
NET VELOCITY 
SUMS 


Yes SUBROUTINE 
NETVQD 


COMPUTE 
NET FLOW 
SUMS 


SUBROUTINE NETVQD 
(ENTRY NETFLO) 


STORE 
INSTANTANEOUS HYDRODYNAMICS 
IF AT PROPER TIME 
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MEAN VELOCITIES 
FOR COMPUTATION 
OF DISPERSION 
OEFFICIENTS 


SUBROUTINE NETVQO 
(ENTRY UVDXDY) 


SUBROUTINE PRINTO 
(ENTRY SAVEQH) 


$ 


DOES 
TIME = TPER 
J] TIDAL PERION 


3 


PUNCH 
HYDRODYNAMICS 


PUNCH H VALUES 
AT SELECTED CELLS 
FOR USE IN FINE 
GRID MODEL 
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IS 
TIME > TMAX 
(TOTAL SIMULATION 
TIME) 


STORE FINAL VALUES 
OF SELECTED FLOWS & 

TIDES FOR DETERMINING 
FINE GRID INPUTS 


SUBROUTINE RITAP 


j ( 


TRANSFORM COARSE 
GRID FLOWS INTO 
FINE GRID EXCIT- 
ING FLOWS 


PRINT HYDRODYNAMICS 
FOR 10 REMAINING 
SPECIFIED CELLS 


SUBROUTINE PRINTO 
(ENTRY PRT11) 


SAVE 
ENDING VALUES OF HYDROS 
DYNAMICS 


SUBROUTINE PRINTO 
(ENTRY SAVEQH) 


COMPUTE 
AND STORE NET 
ELOCITIES 


SUBROUTINE NETVQD 
(ENTRY PVLDEP) 


COMPUTE 
AND STORE NET 
FLOWS 


SUBROUTINE NETVQP 
(ENTRY PFLDEP) 


COMPUTE 
AND STORE DISPERSION 
OEFFICIENTS 


SUBROUTINE NETVQD 
(ENTRY UVDXDY) 


SUBROUTINE PLOTHS 
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FORTRAN LISTING OF PROGRAM HYDTID 


Following is a complete listing of the HYDTID program including 
all of its various subroutines. This listing corresponds to that 
used for simulation of the modified inlet conditions with the jetty 
in place (12 September 1969 and June 1967 simulations). To apply 
the program to the pre-project iniet conditions without the jetty 
(November 1964 simulation), there are four changes which have to be 
made. The sets of cards which have to be interchanged are identified 
in the program listing as Insert 1, Insert 2, Insert 3, and Subroutine 
RITAP. The replacement cards for each of the inserts and a listing 
of Subroutine RITAP, Version I, which applies to the modified inlet 


conditions, follows the complete HYDTID listing. 
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OUTPUT FROM PROGRAM HYDTID 


Two complete sets of sample output from HYDT{D are included 
in this section. Both are for simulation of the verification 
period 0400 to 2100 on 12 September 1969, however the first is 
output from the coarse grid model and the second is from the 
fine grid sub-model. Each set includes exactly the same types 


of information but for different models. 


118 


Coarse Grid Model Output 
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